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Abstract. We study the linear filtering problem for systems driven 
by continuous Gaussian processes Vi and V2 with memory described 
by two parameters. The processes Vj have the virtue that they pos- 
sess stationary increments and simple semimartingale representations 
simultaneously. It allows for straightforward parameter estimations. 
After giving the semimartingale representations of Vj by innovation 
theory, we derive Kalman-Bucy-type filtering equations for the sys- 
tems. We apply the result to the optimal portfolio problem for an 
investor with partial observations. We illustrate the tractability of 
the filtering algorithm by numerical implementations. 



1. Introduction 

In this paper, we use the following Gaussian process (^(t))teR with 
stationary increments as the driving noise process: 

(1.1) V{t) = W{t) - (^J pe-('?+P)(^-")rfW^(M)^ ds {t G R), 

where p and q are real constants such that 

(1.2) < g < 00, -g < p < 00, 

and (W(t))t£-R is a one- dimensional Brownian motion satisfying W{0) = 0. 
The parameters p and q describe the memory of V"(-). In the simplest case 
p = 0, V{-) is reduced to the Brownian motion, i.e., V^(-) = W{-). 

It should be noticed that ()1.1|) is not a semimartingale representation 
of {V(t))o<t<T with respect to the natural filtration J^^{-) of iV{t))o<t<T 
since {W{t))o<t<T is not jF^(-)-adapted. Using the innovation theory as 
described in Liptser and Shiryayev [2Tj and a result in Anh et al. 0, we 
show f Theorem 12. ip that there exists a one- dimensional Brownian motion 
{B(t))o<t<T, called the innovation process, satisfying 

a{B{s) : < s < t) = (t{V{s) : < s < t) (0 < t < T), 
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where l{t, s) is a Volterra kernel given by 





With respect to the natural filtration J^^{-) of B{-), which is equal to 
()1.3|1 gives the semimartingale representation of V^(-). Thus the 
process V{-) has the virtue that it possesses the property of a stationary 
increment process with memory and a simple semimartingale representa- 
tion simultaneously. We know no other process with this kind of proper- 
ties. The two properties of V^(-) become a great advantage, for example, 
in its parameter estimation. 

In jT] and Anh et al. |21IS], the process V^(-) is used as the driving pro- 
cess for a financial market model with memory. The empirical study for 
S&P 500 data in Anh et al. |4J shows that the model captures very well 
the memory effect when the market is stable. The work in these references 
suggests that the process V^(-) can serve as an alternative to Brownian mo- 
tion when the random disturbance exhibits dependence between different 
observations. 

In this paper, we are concerned with the filtering problem of the two- 
dimensional partially observable process {X{t),Y{t))Q<t<T governed by 
the following linear system of equations: 



Here X{-) and Y{-) represent the state and the observation respectively. 
For j = 1,2, the noise process Vj{-) is described by p.ll) with {p,q) = 
{Pj,qj) and W{-) = Wj{-). We assume that the Brownian motions Wi{-) 
and W2{-), whence Vi(-) and V2(-), are independent. The coefficients 
9,a, fi eH with 7^ are known constants, and Xq is a centered Gaussian 
random variable independent of (Vi, V2). 

The basic filtering problem for the linear model ()1.5|) with memory is 
to calculate the conditional expectation E[X(t)|jF^(t)], called the (opti- 
mal) filter, where J^^(-) is the natural filtration of the observation process 
Y{-). In the classical Kalman-Bucy theory (see Kalman |T2|, Kalman and 
Bucy jlSI, Bucy and Joseph Davis and |2I]), Brownian motion is 
used as the driving noise. Attempts have been made to resolve the filter- 
ing problem of dynamical systems with memory by replacing Brownian 
motion by other processes. In Kleptsyna et al. [THl [13 HHj and others, 
fractional Brownian motion was used. Notice that fractional Brownian 



(1.5) 
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motion is not a semimartingale. In the discrete-time setting, autoregres- 
sive processes are used as driving noise (see, e.g., [12], [S] and Jazwinski 
[TT]). Our model may be regarded as a continuous-time analogue of the 
latter since it is shown in pP that V{-) is governed by a continuous-time 
AR(cx))-type equation. 

The Kalman-Bucy filter is a computationally tractable scheme for the 
optimal filter of a Markovian system. We aim to derive a similar effec- 
tive filtering algorithm for the system (jl.5|) which has memory. However, 
rather than considering (jl.5|) itself, we start with a general continuous 
Gaussian process {X(t))o<t<T as the state process and Y defined by 



as the observation process, where is a deterministic function and V{-) 
is a process which is independent of X(-) and given by Using (jl.3|) 

and (jl.4|) . we derive explicit Volterra integral equations for the optimal 
filter f Theorem 13. In the special case ()1.5|) . the integral equations are 
reduced to differential equations, which gives an extension to Kalman- 
Bucy filtering equations f Theorem 13. 4|) . Due to the non-Markovianness of 
the formulation ()1.5|) . it is expected that the resulting filtering equations 
would appear in the integral equation form (cf. Kleptsyna et al. JH]). The 
fact that good Kalman-Bucy-type differential equations can be obtained 
here is due the special properties of (|1.5p . This interesting result does not 
seem to hold for any other formulation where memory is inherent. 

We apply the results to an optimal portfolio problem in a partially 
observable financial market model. More precisely, we consider a stock 
price model that is driven by the solution V{-) to Assuming that 

the investor can observe the stock price but not the drift process, we 
discuss the portfolio optimization problem of maximizing the expected 
logarithmic utility from terminal wealth. To solve this problem, we make 
use of our results on filtering to reduce the problem to the case where 
the drift process is adapted to the observation process. We then use the 
martingale methods (cf. Karatzas and Shreve to work out the explicit 
formula for the optimal portfolio f Theorem 14.11) . 

This paper is organized as follows. In Section El we prove the semi- 
martingale representation ()1.3|) with ()1.4j) for V{-). Section El is the main 
body of this paper. We derive closed form equations for the optimal filter. 
In Section HI we apply the results to finance. In Section we illustrate 
the advantage of V{-) in parameter estimation. Some numerical results on 
Monte Carlo simulation are presented. Finally, in Section ^ we numeri- 
cally compare the performance of our filter with the Kalman-Bucy filter 
in the presence of memory effect. 
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2. Driving noise process with memory 



Let T G (0,00), and let (fi, JF, P) be a complete probability space. For 
a process {A{t))o<t<T, we denote by jF^(t) the P-augmentation of the 
filtration a{A{s) : < s < t), < t < T. 

Let V{-) be the process given by The process (y(t))tgR, is a 

continuous Gaussian process with stationary increments. The aim of this 
section is to prove ()1.3p with ()1.4|) . 

By |211 Theorem 7.16], there exists a one-dimensional Brownian motion 
{B{t))o<t<T, called the innovation process, such that J-'^{t) = J-'^{t), 
< t < T and that 



(2.1) 



V{t) = B{t) 



a{s)ds, (0 < t < T), 



pe 



-{(l+'P){t-u 



u 



(0 < t < T). 



Thus, {V(t))Q^t<T is a jF^-semimartingale with representation (|2.1|) . 

To find a good representation of a(-), we recall the following result from 
121 Example 5.3]: 



with 



k{t, s) = p{2q + p) 



k{t,s)dV{s) (0<t<r) 
(2g + p)e'?" - pe"''" 



(2g + p)2e9* - p2e-9t 



(0 < s < t). 



From the theory of Volterra integral equations, there exists a function 
l{t,s) G L^[0,T]^, called the resolvent of k{t,s), such that, for almost 
every < s < t < T, 



(2.2) 



(see [6, Chapter 4, Section 3] and [TUl Chapter 9]). Using l(t,s), we have 
the following representation of a in terms of the innovation process B: 



l{t, s) — k{t, s) + J l{t,u)k{u, s)du = 0, 
l{t, s) — k{t, s) + / k{t,u)l{u, s)du = 



(2.3) 



a{t) 



l{t,s)dB{s) {0<t<T). 



We shall solve (j2.2p explicitly to obtain l(t,s). 

Theorem 2.1. The expression l\1.4l holds. 
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Proof. We have k{t, s) = a{t)b{s) for < s < t, where, for t G [0,T], 

Fix s G [0,T] and define x{t) = Xgit) and A = by 

x(t) = I b{u)l{u, s)du {s<t<T), A = b{s). 

J s 

Then, from ()2.2|) we obtain 

— {t) + a{t)b{t)x{t) = \a{t)b{t), x{s) = 0. 
The solution x is given by 

x{t)=X- Ae-i'>(")^(")'^", 
whence l{t, s) is obtained as 

l{t, s) = a(t)6(s)e-^X")^(")'^" = k{t, s)e-/X"''')'^« (0 < s < t). 
We have 

K[U, U) = p — 



(2g + p)2e29« 

By the change of variable x{u) = {2q + p)^e^'^" — p^, we obtain 

^2A ri--^i^-io.^^^)^^^^)+^^^ 



so that 



Thus 



xis) X+p^j X{s){x{t) + p"^) 

^ ^ (2g + p)2e2'?« - p2 



(2g + p)2e29^ 



(2g + p)2e2'?« -p2 
or p.4|) . as desired. 



3. Filtering equations 



3.1. General result. Let (X(t), ?7(t))o<t<T be a two-dimensional cen- 
tered continuous Gaussian process. The process X(-) represents the state 
process, while [/(■) is another process which is related to the dynamics of 
X(-). 

Let (i?(t))o<t<T be a one-dimensional Brownian motion that is inde- 
pendent of {X,U). We define the processes V{-) and a{-) by ()2.1|) and 
(j2.3|) . respectively. In this subsection, we assume that l(t, s) in (j2.3j) is an 
arbitrary Volterra-type bounded measurable function (i.e., l{t,s) = for 
s > t) rather than the special form (jl.4p . Of course, the function l{t, s) in 
p.4j) satisfies this assumption. We consider the observation process Y{-) 
satisfying 

Y{t)= [ fi{s)X{s)ds + V{t) (0<t<T), 







where /x(-) is a bounded measurable deterministic function on [0,T] such 
that /i(t) ^ for < t < T. 

As in Sectional let JF^ be the augmented filtration generated by Y{-). 
For d-dimensional column vector processes (A(t))o<t<T and (C(t))o<i<T, 
we write 

A{t) =E[A{t)\J^''{t)] (0<t<T), 
TAcit, s) = E[A{t)C*{s)] (0 < s < t < T), 

where * denotes the transposition. Notice that TAcit, s) G R'^^'^. 
We put 

Z{t) = iX{t),U{t),a{t)r (0<t<T), 
and define the error matrix P{t,s) G R^^^ by 

P{t, s) = E[Z{t){Z{s) - Z{s)y] (0 < s < t < T). 

The next theorem gives an answer to the filtering problem for the par- 
tially observable process {X{t),Y{t))o<t<T- It turns out that this will be 
useful in the filtering problem for ()1.5|) for example. 

Theorem 3.1. The filter Z(-) satisfies the stochastic integral equation 

(3.1) Z{t)= [ {P{t,s)+D{t,s)}a{s){dY{s)-a*{s)Z{s)ds}, 
Jo 

and the error matrix P{t, s) is the unique solution to the following matrix 
Riccati-type integral equation such that P{t, t) is symmetric and nonneg- 
ative definite for <t <T: 



P{t,s) = - {Pit,r)+Dit,r)}air)a*ir){Pis,r) + Dis,r)rdr 
[6.Z) Jo 

+ Tzzit,s) (0<s<t<T), 
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where 



0' 

D{t,s) = { 

At,s)/fi{s) 0, 



as 



/i(s) 


-1 



Proof. Since (X, U) is independent of B, {X, U, a, Y) forms a Gaussian 
system. We have 

ft 



Y{t) 



{fi{s)X{s) - a{s)}ds + B{t). 



Thus we can define the innovation process /(■) by 

/(t) = Y{t) - [ {fi{s)X{s)ds - a{s)) (0 < t < T), 
Jo 

which is a Brownian motion satisfying JF^ = JF^ (cf. |2Tl Theorem 7.16]). 
Notice that /(■) can be written as 



(3.3) 



lit) 



{Z{s) - Z{s)ya{s)ds + B{t). 



By [21, Corollary of Theorem 7.16], there exists an R^- valued Volterra- 
type function F{t,s) = {Fi{t, s), F2{t, s), Fs^t, s))* on [0,T]2 such that 



[ \F{t, s)\'^ds < +00 (0<t<T), 
Jo 



(3.4) 



Z{t) 



F{t,s)dl{s) (0<t<T), 



where | ■ | denotes the Euclidean norm 

Now let g(t) = {giif) , g2{t) , g^{t)) be an arbitrary bounded measurable 
row-vector function on [0,T]. Then, for t G [0,T], 

E [ / g{s)dl{s) ■ {Z{t) - Z{t)) = 0. 



From this, ()3.3|) . ()3.4|) and the fact that {X,U) and B are independent, 
we have 



[ g{s)F{t,s)ds = E [ g{s)dl{s) ■ Z{t) 
Jo Jo 



E 



g{s){{Z{s) - Z{s)ya{s)ds + dB{s)} ■ Z{t) 



g{s)¥.[Z{t){Z{s)-Z{s)Y]a{s)ds+ / g^{s)l{t, s)ds 



g{s)P{t, s)a{s)ds + / g{s)D{t, s)a{s)ds. 

) 
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Since g{-) is arbitrary, we deduce that F(t, s) = {P{t, s) + D{t, s))a{s) or 

(3.5) Z{t)= [ {P{t,s) + D{t,s))a{s)dl{s) (0 < t < T). 
Jo 

The SDE ()3.H) follows from ()3.5|) and the representation 

I{t) = Y{t) - [ a*{s)Z{s)ds {0<t< T). 
Jo 

The equation ()3.2|) follows from ()3.5|) and the equality 

Pit, s) = E[Z{t)Z*{s)] - E[Z{t)Z*{s)]. 

The matrix P{t,t) is clearly symmetric and nonnegative definite. Finally, 
the uniqueness of the solution to ()3.2j) follows from Proposition 13 .21 below. 

□ 

Proposition 3.2. The solution P{t, s) to the matrix integral equation 
1)5'. such that P{t, t) is symmetric and nonnegative definite for <t <T 
is unique. 

Proof. By coninuity, there exists a positive constant C{T) such that 
\\Tzz{t,s)\\ <C{T) (0<s<t<T), 

where ||v4|| := (trace(v4M))^/^ for A e R^""^. Let P be a solution to (IT^ 
such that P{t,t) is symmetric and nonnegative definite for t G [0,T]. We 
put Q{t, s) = P{t, s) + D{t, s). Then ((23) with s = t is 

rzz{t,t)-P{t,t)= [ Q{t,r)a{r)a*{r)Q*{t,r)dr. 



t pt 

2, 







From this, we have 

\Q{t,r)a{r)\'^dr = I tr ace {Q{t,r)a{r) a* (r) Q* {t, r)} dr 
Jo 

= tiaceiT zz{t,t) - P{t,t)) < tTace{Tzz{t,t)) 
<V3\\Tzz{t,t)\\<V3C{T). 
Therefore, ||P(t, is at most 

\\rzz{t,s)\\ + [ \\Q{t,r)a{r)a*{r)Q*{s,r)\\dr 



<\\Tzz{t,s)\\ + (^j\Q{t,r)a{r)\'dry ■ (^j^ \Q{s,r)a{r)\'dr^ ' 
< (1 + y3)C(T). 

Let Pi and P2 be two solutions of ()3.2|1 . We define s) = Piit, s) + 
D(t,s) for i = 1,2. We put Pi(t,s) = for s > t and i = 1,2. Since /i 



and / are bounded, it follows from the above estimate that there exists a 
positive constant K{T) satisfying 

\\a{s)a*{s)Qi{t,s)\\ < K{T) iO<s<t<T, i = 1,2). 

It follows that 

\\Qi{t, r)a{r)a*{r)Qi{s, r)* - Q^it, r)a{r)a\r)Q2{s, r)*|| 

< ||gi(t, r)a{r)a\r){Q,{s, r)* - Q^is, r)*) || 

+ II iQiit, r) - Q^it, r))air)a*ir)Q,is, r)*|| 

< 2K{T)\\Qr{s,r)-Q2{t,r)\\ = 2K{T)\\P,{s,r) - P2{t,r)\\. 
From this and (j3.2j) . we obtain 

sup \\P,{t,s)-P2{t,s)\\<2K{T) [ sup \\P,{t,r)-P2{t,r)\\dr. 

0<t<T Jo 0<t<T 

Therefore, Gronwall's lemma gives 

sup ||Pi(t,s) -P2(t,s)|| = (0<s<T). 

0<t<T 

Thus the uniqueness follows. □ 

Remark 1. We consider the case in which a = and the state process 
X{-) is the Ornstein-Uhlenbeck process satisfying 

dX{t) = eX{t)dt + adW{t), X{0) = 0, 

where 0, a 7^ and W{-) is a one-dimensional Brownian motion that is 
independent of -B(-). We also assume that /i(-) = fi, a constant. Then 

X{t) = a e^(*-")rfiy(M) and 

rxx(t, s) = ^(e'^*+^^ - e^(*-^)) (0 < . < t < T). 
By Theorem 13. If we have 

(3.6) X{t)= [ fiPxx{t,s){dY{s) - fiX{s)ds), 

Jo 

(3.7) Pxxit,s) = Txxit,s)- [ fi^Pxxit,r)Pxx{s,r)dr, 



where Pxx{t,s) = E[X{t){X{s) - X{s))] for < s < t < T. Let 
jF(t), < t < T, be the P-augmentation of the filtration generated by 
{W{t),B{t))o<t<T. Then Pxx{t,s) is 

E[X{t){X{s) - X{s))] = E[E[X{t)\J^{s)]{X{s) - X{s))] 

= E[e'^'-'^X{s){X{s) - X{s))] = e'^'-'h{s) 

with 7(s) = E[X{s){X{s) - X{s))]. Thus (jHSl) is reduced to 

(3.8) dX{t) = {6- ^'^^{t))X{t)dt + /i7(t)rfy(t). 



Differentiating (|3.7p in s = t, we get 

(3.9) ^ = a^ + 2^7(t)-A(tf. 

The equations ()3.8|) and ()3.9|) are the well-known Kalman-Bucy filtering 
equations (see [ISl, 0, 0, CH and EH)- 

3.2. Linear systems with memory. Let (X(t), F(t))o<i<T be the two- 
dimensional partially observable process satisfying (|1.5p . For j = 1,2, 
the noise process Vj{-) there is described by with W{-) = Wj{-) 

and {p,q) = ipjAj) satisfying (jl.2p . The Brownian motions Wi{-) and 
W2{-)^ whence Vi(-) and V2(-), are independent. In ()1.5|) . the coefficients 
^, 0", /i G R with yU 7^ are known constants and the initial value Xq is a 
Gaussian random variable that is independent of (Vi, V2). The processes 
X(-) and F(-) represent the state and the observation, respectively. 
By Theorem 12.11 we have Vj{t) = Bj{t) — aj{s)ds with aj{t) = 

lj{t, s)dBj(s) for j = 1,2 and < t < T, where -Bi(-) and B2{-) are two 
independent Brownian motions, and, for < s < t < T, 

We put lj{t) = lj{t, t) for j = 1, 2, that is. 

It holds that /^(t, s) = e~^^^''^~^Hj{s), where rj = pj + qj. 

We denote by J^{t), < t < T, the P-augmentation of the filtration 

CT{Xo,{V,{s),V2{s))o<s<t),0<t<T. 

Lemma 3.3. For < s < t < T , we have 

E[X{t)\J^{s)] = e^^'-'^X{s) - ab{t - s)ai(s), 

where 

'(e^*-e-'-i*)/(^ + ri) (^ + ri^O), 



m 



n = 0). 



t 

h{t-u)li{u)dBi{u). 



Proof. For t G [0,r], j^^e^^^-'^ai{s)ds is 

Since = e^*Xo + (x /g e^(*-")rfV^i(M) or 

X(t) = + a / e^(*-")rf5i(M) - a /" e^(*-")ai(u)rfu, 

Jo Jo 



E[X{t)\J^{s)] with s<tis equal to 

e^'Xo + ae^^'-'^ [ e^^'-''^dBi{u) - a [ b{t - u)h{u)dBi{u) 
Jo Jo 

= e^(*-^)X(s) -a r {b{t -u)- e'^'-'\s - u)] h{u)dB,{u). 





However, by elementary calculation, we have 

b{t -u)- e^(*~")6(s -u) = b{t - s)e-'''^'-''^ {0<u<s<t). 
Thus the lemma follows. □ 
We put, for < t < T, 

/ 0\ / fi 

F = \ n 1, D{t) = , a = 
ra/ \f^-%{t) 0/ 

ah{t) 0\ f -e a 

G{t) = I cxk{t) Ijit) , H{t) = ri 

0/ \/i/2(t) r2-/2(t), 

We also put 

Z{t) = {X{t),ai{t),a2{t)y {0<t< T). 
Recall that rzz(O) = E[Z(0)Z*(0)] and that 

Z{t) =¥.[Z{t)\J^^{t)] {0<t<T). 
We define the error matrix P{t) G R^^'^ by 

P(t) = E[Z{t){Z{t) - Z{t))*] {0<t< T). 
Here is the solution to the optimal filtering problem for (jl.Sp . 
Theorem 3.4. The filter Z{-) satisfies the stochastic differential equation 

dZ(t) = -{F + (P(t) + D{t))aa*}Z{t)dt 

+ {P{t) + D{t))adY{t) (0<t<T), 

with Z{0) = (E[Xo], 0, 0)*, and P(-) follows the matrix Riccati equation 

(3 12) ^ = ^(^) - H(t)P{t) - Pmity - P{t)aa*P{t) 

(0 < t < T) 

wzth P,,{0) = 5a5,iE[(Xo)2] /or 2, J = 1,2,3. 
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Proof. For < s < t < T, we put 

P{t,s)=E[Zit)iZis)-Zis)y]. 
Then we have P(t) = P(t, t). We also put, for < s < t < T, 
D{t,s) = e-"2(*-^)D(s), 

Q{t, s) = P{t, s) + D{t, s), Q{s) = Q{s, s) = P{s) + D{s). 

By 

P{t, s) = E[E[Z(t)|^(s)](Z(s) - Zis))*] 
and Lemma EISl we have P{t, s) = M{t — s)P{s), where 

/e^* -ab{t) 
M{t) = e-'i* 

\ e-^2t^ 

with as in Lemma We also see that D{t, s) = M{t — s)D{s). Com- 
bining, Q{t, s) = M{t - s)Q{s). However, M{t) = e'^^ since dM{t)/dt = 
—FM{t) and M(0) is the unit matrix. Thus we obtain 

(3.13) Q(t,s) = e-(*-^)^Q(s). 

From ()3.13|) and Theorem 13 . II with U = ai and a = a2, it follows that 

(3.14) Z{t) = [ e-^'-'^^Q{s)a{dY{s) - a*Z{s)ds}, 

Jo 

(3.15) P(t) = Tzz{t)- f e-(*-")^g(M)aa*Q*(u)e-(*-")^*ciM. 

The SDE flTTni follows easily from (jSHD. 
Differentiating both sides of ()3.15|) yields 

P = tzz + F{Tzz -P) + {Tzz - P)F* - Qaa*Q* 

= tzz + FTzz + TzzF* - Daa*D - HP - PH* - Paa*P. 

Since 

dZ{t) = -FZ{t)dt + dR{t) 

with 

R{t) = {oBiit), k{s)dBi{s), l2{s)dB2{s) 
we see by integration by parts that Z{t)Z{t)* — Z{0)Z{0)* is equal to 

Z{s)dZ{s)*+ [ dZ{s)Z{s)* + E[R{t)R*{t)]. 
Jo 
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It follows that T zzif) — rzz(O) is equal to 



E 



Z{s)dZ{s)* 



E 



dZ{s)Z{sy 



+ E[i?(t)i?*(t)] 



( Tzzis)F*ds- [ FTzzis)ds + E[R{t)R*{t)]. 
Jo 



Tzz{t)+ / FTzz{s)ds+ / Tzz{s)F''ds-Vzz{^) 



Thus 



= E[R{t)R*{t)] 
This and (pnH|) yield (PTT^ . 



+ L'(s)aa*L)*(s))rfs. 



□ 



Remark 2. We can easily extend Theorem IH.4I to a more general setting 
where the functions s) take the form s) = e^^*~''''(7(s). 

Corollary 3.5. We assume that p2 = 0, i.e., V2(-) = W2{-). Let Z{t) = 
{X{t),ai{t)y and P{t) = E[Z{t){Z{t) - e R^^^^ Then the filter 

Z{-) and the error matrix P(-) satisfy, respectively, 

dZ{t) = -{F + P{t)aa*}Z{t)dt + P{t)adY{t), Z{0) = (E[Xo], 0)*, 



dPjt) 

dt 
where 



G{t) - FP{t) - P{t)F* - P{t)aa*P{t), P,,(0) = 5a5,MXo)% 

ah{ty 



-9 a 
rij ' 



, Git) 



ah{t) h{tf 



Corollary 3.6. We assume that pi = 0, i.e., Vi{-) = Wi{-). Let Z{t) = 
(X(t),a2(t))* and P{t) = E[Z{t){Z{t) - Z{t)y] E R^^^^ Then the filter 
Z{-) and the error matrix P(-) satisfy, respectively, 

dZ{t) = -{F + {P{t) + D{tyaa*}Z{t)dt + (P(t) + D{tyadY{t), 
dP{t) 



dt 



G - H{t)P{t) - P{t)H{ty - P{t)aa*P{t), 



with Z(0) = (E[Xo],0)* and Py(0) = (5,i(5jiE[(Xo)2], wh 



ere 



G 



-e 

r2 







, D{t) 











-e 

fihit) r2 - kit) J ■ 
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Example 3.7. We consider the estimation problem of the value of a signal 
p from the observation process (l^(t))o<t<T given by 



dY{t) = pdt + dV{t), F(0) = 0, 

where V{-) and «(■) are as in Section |21 We assume that p is a random 
variable having the normal distribution with mean zero and variance v^. 
This is the special case = cr = of the setting of Corollary 13.61 Let 
r = p + q and /(■) be as above. Let H{t) and a be as in Corollarv 13.61 
with p = 1 and ^ = 0. We define P(t) = (Pii(t))i<i,i<2 by P{t) = 
E[Z*{t){Z{t) - Z{t))] with Z{t) = ip,a{t)y. Then, by Corollarv ITHl the 
filter (p(t),a(t)) satisfies 

dp{t) = -{Pn{t) - Pi2{t)}{p{t) - a{t)}dt + {Pn{t) - Pu{t)}dY{t), 
da{t) = [{P2i(t) - P22{t) + m - r}a{t) - {P^iit) - P22(t) + /(t)}p(t)]rft 

+ {P2l{t)-P22{t) + l{t)}dY{t) 

with (p(0),a(0)) = (E[Xo],0), and the error matrix P(-) follows 

By the linearization method as described in Chapter 5], we can solve 
the equation for P(-) to obtain 

( 1 -mim 



^^^^ 1 + v^r]{t) + v^m(l>it) V-<^(^)/^W v^<Pit?l^it? 
where 

i){t) = exp ( / (r - l{s))ds\ , = /" l{s)ip{s)ds 







The analytical forms of ijj, 0, ^ and r/ can be derived. We omit the details. 

4. Application to finance 

In this section, we apply the results in the previous section to the prob- 
lem of expected utility maximization for an investor with partial observa- 
tions. 

Let {Vj(t))o<t<T, {cijit))o<t<T, j = 1, 2, be as in Section|21 In particular, 
Vi{-) and V2(-) are independent. We consider the financial market consist- 
ing of a share of the money market with price S'o(t) at time t G [0, T] and 
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a stock with price Sit) at time t G [0, T]. The stock price S{-) is governed 
by the stochastic differential equation 

(4.1) dS{t) = S{t){U{t)dt + vdV2{t)}, S{0) = sq, 

where Sq and r] are positive constants and U{-) is a Gaussian process 
following 

(4.2) dU{t) = {6 + eU{t)}dt + adViit), U{0) = p. 

The parameters 6, 6 and a are constants, and p is a Gaussian random 
variable that is independent of (Vi, V2). For simplicity, we assume that 

5o(-) = l, r/=l, 5 = 0. 

Let J-'{t), < t < T, be the P- augment at ion of the filtration generated 
by (Vi(s), V2(s))o<s<t and p. Then f/(-) is jF-adapted but not jF-^-adapted; 
recall from Section|21that JF"^ is the augmented filtration generated by S{-). 
Suppose that we can observe neither the disturbance process V2(-) nor the 
drift process f/(-) but only the price process S{-). Thus only ^'^-adapted 
processes are observable. 

In many references such as f?^, fH^, and jH], the partially observable 
model described by ()4.1|1 and ()4.2|1 with V^'s replaced by Brownian mo- 
tions, i.e., Vj = Bj, is studied. 

We consider the following expected logarithmic utility maximization 
from terminal wealth: for given initial capital x G (0, 00), 

(4.3) maximize E[log(X^'^(T))] over all it G A{x), 
where 

7r(-) is R-valued, jF"^-progressively measurable,! 
• j^'K\t)dt < 00, X"'-(t) > (0 < t < T) a.s. J ' 

and 

(4.4) X--{t)=x+ f^,dS{u). 

Jo S{u) 

The value 7i{t) is the dollar amount invested in the stock at time t, whence 
7i{t)/S{t) is the number of units of stock held at time t. The process 
X^''^{-) is the wealth process associated with the self- financing portfolio 
determined uniquely by 7r(-). 

An analogue of the problem ()4.3|) for full observations is solved in |2]. 
For related work, see jT3], ||20j and the references therein. We solve 
the problem ()4.3|) by combining the results above on filtering and the 
martingale method as described in |14j . 

Consider the process (l^(t))o<t<T defined by 

Y{t)= [ U{s)ds + V2{t) =B2{t)+ [ {U{s) - a2{s))ds, 
Jo Jo 
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which we regard as the observation process. Let (/(t))o<t<T be the inno- 
vation process associated with Y{-) that is given by 

I{t) = Y{t) - [ {U{s) - a2{s))ds, 
Jo 

where A{t) = E[y4(t)|jF^(t)] as in the previous sections. The innovation 
process /(■) is a JF^-Brownian motion satisfying J-'^ = JF^ = J^^ . 
Let (L(t))o<t<T be the exponential jF-martingale given by 

m = exp i^-j\u{s) - a2{s))dB,{s) - ^ j\u{s) - «2(s))'rfs| 

= exp |- ^ {U{s) - a2{s))dY{s) + ^ ^ {U{s) - a2{s)fds 
By EH Chapter 7], we find that, for t G [0,T], 

Lit) = exp |- ^ (f/(s) - a2{s))dY{s) + ^ ^ (f7(s) - a2{s)fds 

= exp |- ^ {U{s) - a2{s))dl{s) - ^ ^ {U{s) - a2{s)fds 

The process (L(t))o<t<T is a jF^-martingale. 

For X G (0, cxd) and tt G A{x), we see from Ito formula that the pro- 
cess {L{t)X^'^{t)) is a local jF^-martingale, whence a jF^-supermartingale 
since X^''^(-) is nonnegative. It follows that, for x G (0, oo), vr G ^(x), 
and y G (0, oo), 

(4.5) E[log(X^'-(T))] < E[log(X^'-(r)) - |/L(r)X^''-(T)] + yx 

<E[log{l/{yL{T))}-l]+yx, 
where, in the second inequality, we have used 

^og{z) -yz < log(l/y) - 1 (y, z G (0, oo)). 
The equalities in ()4.5p hold if and only if 

(4.6) X^'"(T) = x/L{T) a.s. 

Thus the portfolio process 7r(-) satisfying (j4.6j) is optimal. 
Put 

7ro(t) = xiUit) - a2{t))/L{t) {0 < t < T). 
Since (x/L)(0) = x and 
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we see from (j4.4|) that the process vro(-) satisfies (j4.6|) . whence it is the 
desired optimal optimal portfolio process. It also satisfies 

^L = f/(t)-a,(t) (o<t<r). 

We put 

a = (1,0,-1)*, Z{t) = {U{t),ai{t),a2{t)r {t e [0,T]). 

We define the error matrix P{t) e R^^^ E[Z{t){Z{t) - Z{t))*]. Com- 
bining the results above with Theorem 13.41 which describes the dynamics 
of U and 02, we obtain the next theorem. 

Theorem 4.1. The optimal portfolio process ttq for the problem l\4-^ is 
given by 

7ro(t) = xa*Z{t)/L{t) (0 < t < T) 

and satisfies 

X^'-(T) = a:/L(T), = (0 < t < T). 

The filter Z(-) follows 

dZ{t) = -{F + (P(t) + D{t))aa*}Z{t)dt + (P(t) + D{t))adY{t) 

with Z{0) = (E[p],0,0)*, and the error matrix P(-) satisfies the matrix 
Riccati equation 

dP{t) 



dt 



G{t) - H{t)P{t) - P{t)H{t)* - P{t)aa*P{t) 



with Pij{0) = 5ii5jiE[p2] = 1,2,3), where F, D{t), G{t) andH{t) are 
as in Theorem with fi = 1. 

5. Parameter estimation 

Let V^(-) be the process given by (jl.lj) . We can estimate the values of 
the parameters p and q there from given data of V{-) by a least-squares 
approach (jlj). In fact, since V^(-) is a stationary increment process, the 
variance of V(t) — V{s) is a function in t — s. To be precise, 

^ VaT{V{t) -V{s)) = U{t~ s) {0<s<t), 



t-s 

where the function U{t) = U (t; p, q) is given by 

[p + qy [p + qy^ t 
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Suppose that for tj = jAt, j = 1,. . . ,N, the value of V(tj) is Vj. For 
simphcity, we assume that At = 1. The unbiased estimation Uj that 
corresponds to U(tj) is given by 

N-j 

where rrij is the mean of Vi+j — fj's: 

^ N-j 

1=1 

Fitting {U{tj]p,q)} to {uj} by least squares, we obtain the desired esti- 
mated values of p and q. 

For example, we produce sample values f i, f2, . . . , f looo ioi V{-) with 
(p, q) = (0.5, 0.3) by a Monte Carlo simulation. We use this data to 
numerically calculate the values po and go of P and q, respectively, that 
minimize 

30 

J2iu{tr,p,q)-u,)^. 

i=i 

It turns out that po = 0.5049 and go = 0.2915. In Figure EIH we plot 
{U{tj;p,q)} (theory), {uj} (sample) and {U(tj;po,qo)} (fitted). It is seen 
that the fitted curve follows the theoretical curve reasonably well. 

We extend the approach above to that for the estimation of the param- 
eters p, g, 9 and a in 

dX{t) = -eX{t)dt + adV{t), X(0) = Xo, 

where 9,a E (0,oo), the process V{-) is given by as above, and the 
initial value Xq is independent of (^(t))o<t<T and satisfies IE[-^o] < 
The solution X(-) is the following Ornstein-Uhlenbeck-type process with 
memory: 

(5.2) X{t) = e-'% + re-'(*-")(iy(M) {t G [0,T]). 

Jo 

Put 



^(t) := [ e(^-P-«)"c/M (0 < t < T). 
Jo 



Proposition 5.1. We have 

(5.3) j^E[{X{t)-e-'^'-'^X{s)f]=H{t-s) (0 < s < t < T), 
where, for < t < T , the function H{t) = H{t; p, g, 6, cr) is given by 
H(t) = a^{l- P(2g + p) | 1 - e'^^* ^ a^'^q + p)e-2^V W 



{p + q){0 + p + q) \ 2et {p + q){e + p + q)t ' 
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Figure 5.1. Plotting of the function v{t), the sample data 
and the fitted function vo(t). 



Proof. By ()5.2p . the left-hand side of ()5.3|) is equal to 



t-s 



E 



We put Cu = pe-(P+'')"/(o,oo)(M) for m e R. By Proposition 3.2 in j2], 

e^'^dV{u) is given by 



OO / ft 



.udrj dW{u) 



V"- f e^''cr-udr^dW{u)- 



Thus E[(/^*e^"rfV(M))2] is equal to 



s / rt 



e'^^'cr-udrj dW{u). 
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e^^c{r — M)(ir I du. 



By integration by parts and the equalities 

J u 

we obtain the desired result. □ 

Suppose that for tj = jAt, j = 1, . . . ,N, the value of X{tj) is Xj. We 
assume At = 1 for simplicity. The estimation hj{d) that corresponds to 
H{tj;p,q,9,a) is given by 

^ N~j 

where mj{9) is the mean of Xj+j — e^^^Xi, i = 1, . . . , N — j: 

^ N-j 

i=i 

Fitting {H(tj;p,q,9,a) — hj{9)} to {0} by least squares, we obtain the 
desired estimated values of p, g, 9 and a. 

For example, we produce sample values Xi, X2, . . . , Xiooo 

for X(-) with 

a) = (0.2, 1.5, 0.8, 1.0) 

by a Monte Carlo simulation. We use this data to numerically calculate 
the values po; %i 9q and do of p, g, 9 and a, repectively, that minimize 

30 

Y,{H{t,-p,q,9,a)-h,{9)f. 
i=i 

It turns out that 

(po,go,^o,fTo) = (0.1910,1.5382,0.8354,1.0184). 

In Figure we plot {H{tj;p,q,9,a)} (theory), {hj{9o)} (sample with 
estimated 9) and {H(tj;pQ,qo,9o,(7o)} (fitted). It is seen that the fitted 
curve follows closely the theoretical curve. 

6. Simulation 

As we have seen, the process V^(-) described by (jl.ll) has both station- 
ary increments and a simple semimartingale representation as Brownian 
motion does, and it reduces to Brownian motion when p = 0. In this 
sense, we may see V^(-) as a generalized Brownian motion. Since V{-) is 
non-Mar kovian unless p = 0, we have now a wide choice for designing 
models driven by either white or colored noise. 
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5 10 15 20 25 30 

Time lag t 

Figure 5.2. Plotting of the function h{t), the sample data 
h{t) with estimated 6 and the fitted function ho{t). 



In this section, we compare the performance of the optimal filter with 
the Kalman-Bucy filter in the presence of colored noise. We consider 
the partially observable process {X{t),Y(t))o<t<T governed by (jl.Sp with 
Xq = 0. Suppose that an engineer uses the conventional Markovian model 

dX'{t) = eX'{t)dt + adBi{t), X'{0) = 

to describe the non-Markovian system process X{-). Then he will be led 
to use the Kalman-Bucy filter X{-) governed by 

dX{t) = {6- iJ,^-f{t))X{t)dt + iJ.-f{t)dY{t), X(0) = 0, 
where 7(-) is the solution to 

^=a' + 29j{t) - 7(0) = 



(see (|3.8|) and (|3.9|) ). instead of the right optimal filter X as described by 
Theorem 13.41 
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We adopt the following parameters: 

T=10, At = 0.01, 

ti = iAt {i = 1, . 
a = l, e = -2, 



= T/At 
..,N), 
/i = 5. 



1000, 



Let n G {1,2,..., 100}. For the n-th run of Monte Carlo simulation, we 
sample . . . , Xnitj^) for X(-). Let and n = 1, . . . , 100, be 

the Kalman-Bucy filter and the optimal filter respectively. For Un = Xn 
or Xn, we consider the average error norm 



AEN := 

and the average error 



\ 



N 100 



lOOiV 



i=l n=l 



AE{ti 



\ 



^ 100 



Un{ti)y 



1, 



,iV). 



n=l 



In Table 16.11 we show the resulting AEN's of {x^} and for the 
following five sets of 9 = (pi, qi,p2, ^2)'- 

61 = (0.2,0.3,0.5,0.2), 

62 = (5.2,0.3,-0.5,0.6), 

63 = (0.0,1.0,5.8,0.7), 

64 = (5.4,0.8,0.0,1.0), 

65 = (5.1,2.3,4.9,1.3). 

We see that there are clear differences between the two filters in the cases 
02 and 04. We notice that, in these two cases, pi is large than the 
parameters p2 and q2- In Figure lUTTl we compare the graphs of AE(-) for 
the two filters in the case B = 62. It is seen that the error of the optimal 
filter is consistently smaller than that of the Kalman-Bucy filter. 
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Optimal filter 


Kalman-Bucy filter 


ei 


0.5663 


0.5667 


62 


0.4620 


0.5756 


63 


0.5136 


0.5167 


64 


0.4487 


0.5196 


©5 


0.4294 


0.4524 



Table 6.1. Comparison of AEN's 
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Optimal filter 
Kaiman-Bucy filter 




Figure 6.1. Plotting of AE(-) for the optimal and Kalman- 
Bucy filters with noise parameter = 02 
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